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ABSTRACT 


Apparatus  was  built  for  the  purpose  of  measuring 
the  absorption  of  sound  by  the  pulse  technique  in  gases 
at  low  temperatures .  The  fixed  frequency  of  operation  was 
3  Mc/sec,  and  by  varying  the  pressure  of  gas,,  absorption 
measurements  were  made  at  a  f requency-to-pressure  ratio 
ranging  from  3  to  45  Mc/sec-atm.  No  reports  appear  in  the 
literature  about  the  use  of  the  pulse  technique  at  low 
temperatures . 

The  instrument  was  tested  for  proper  operation  by 
measuring  the  absorption  of  sound  in  helium  both  at  room 
temperature  and  at  liquid  nitrogen  temperature. 

The  absorption  of  ultrasonic  waves  was  measured, 
at  liquid  nitrogen  temperature,  in  normal  hydrogen  and  in 
a  para.-enriched  mixture  of  hydrogen. 

The  theory  of  sound  absorption  in  fluids  is 
summarized?  the  first  part  of  the  treatment  deals  with  the 
classical  absorption  due  to  viscosity  and  thermal  conductivity, 
while  the  second  part  deals  with  relaxational  absorption  due 
to  a  slowness  of  energy  exchange  between  the  translational 
and  the  internal  degrees  of  freedom  of  a  gas  molecule.  The 
review  of  the  theory  given  here  is  more  general  than  is 
customarily  given  on  work  in  gases. 

The  results  are  discussed  and  compared  with  similar 
results  obtained  by  different  techniques  at  different 


. 


temperatures . 

The  relaxation  time  for  the  0 - >  2  transition 

in  hydrogen  at  76.7  °K  is  found  to  be  approximately  2.1 
-8 

X  10  sec.  Measurements  both  in  normal  hydrogen  and  in 
para -enri ched  hydrogen  yield  the  same  value  of  the  relaxa¬ 
tion  time;  this  is  in  contrast  to  the  results  of  Sluijter., 
Knaap  and  Beenakker.  The  theories  for  calculation  of 
relaxation  times  are  not  precise  enough  to  shed  light  on 


this  discrepancy. 
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CHAPTER  I 


THEORY 


Introduction 

The  propagation  of  sound  in  fluids  has  been  a 
subject  of  interest  for  well  over  a  century.  The  first 
theoretical  treatment  of  importance  on  the  subject,  dates 
back  to  1845.  It  is  due  to  Stokes,  who  investigated  the 
effects  of  viscosity  on  the  propagation  of  sound.  Later, 
in  1868,  Kirchhoff  showed  that  heat  conduction  in  fluids 
gives  rise  to  an  absorption  of  sound  which  is  of  the  same 
order  of  magnitude  as  that  due  to  viscosity.  A  few  authors, 
among  whom  Lord  Rayleigh  is  prominent,  have  also  evaluated 
the  effect  of  heat  radiation  on  the  propagation  of  sound  in 
fluids;  their  estimates  showed  however,  that  radiation  plays 
little  role  in  gases,  except  perhaps,  at  very  low  pressures 
and  very  high  temperatures. 

Early  experimental  works  consistently  yielded 
results  which  indicated  that  the  absorption  of  sound  in  gases 
other  than  monatomic  gases,  was  much  in  excess  of  that  predicted 
by  the  so-called  classical  theory. 
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In  1928,  Herzfeld  and  Rice  suggested  that  the 
excess  absorption  of  sound  in  fluids  might  be  due  to  a 
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phenomenon  that  they  termed  "thermal  relaxation."  Much 
refined  in  the  following  years ,  this  theory  has  gained  wide 
acceptance,  and  accounts  for  a  great  number  of  experimental 
results.  Sound  absorption  data  in  fluids  can  yield  very 
valuable  information  about  the  nature  of  the  interaction 
between  colliding  particles,  since  the  energy  transfer 
depends  strongly  upon  the  non-spherical  part  of  that  interaction 

A  recent  treatment  of  the  subject  by  Bauer,1  based 
on  the  principles  of  the  thermodynamics  of  irreversible  pro¬ 
cesses,  is  of  particular  interest  in  view  of  its  relatively 
simple,  and  surely  very  elegant  solution  of  such  complicated 
problems  as  multiple  relaxation  processes  in  a  three-level 
gas,  or  in  mixtures  of  two  gases  with  excitable  states. 

Therefore,  a  unified  treatment,  using  the  methods 
of  the  thermodynamics  of  irreversible  processes,  as  developed 
by  Meixner,  will  yield  the  "classical"  sound  absorption 
equations,  also  known  as  the  Stokes -Navier  equations,  as  well 
as  those  for  the  "anomalous"  sound  absorption,  due  to  relaxation 
Meixner  proved  that,  to  a  first  approximation,  the  contributions 
to  the  absorption  of  sound  from  the  transport  phenomena  and 
those  from  the  internal  transitions,  are  additive.  The  two 
problems  can  thus  be  treated  separately. 


lo  a j 

’  ■ 
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Before  dealing  especially  with  the  topic  of  sound 
absorption  in  gases,  a  summary  of  the  fundamentals  of  the 
method  used  is  given. 

1 .Fundamentals  of  the  Phenomenological  Thermodynamics  of 
Irreversible  Processes 

The  usual  inequality  of  ordinary  thermodynamics 
relates  the  increase  in  entropy  dS  of  a  system,  to  its 
absolute  temperature  T,  and  to  the  heat  supplied  to  it,  dQ, 
namely 

dS  >  .  (1-1) 

(l-l)  can  be  made  an  equality  by  writing  as  follows': 

dS  d.S  =  d  S  +  d.S,  (1-2) 

in  which  case  the  differential  deS  represents  the  entropy 
supplied  to  the  system  by  its  surroundings,  and  d^S  represents 
the  entropy  produced  inside  the  system;  d^S  is  always  positive. 

One  of  the  most  important  objectives  of  the  thermo¬ 
dynamics  of  irreversible  phenomena  is  to  relate  d^S  to  the 
various  irreversible  phenomena  which  may  occur  within  the 
system.  In  order  to  do  so,  it  is  convenient  to  use  the  entropy 
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production  per  unit  time  and  unit  volume  -  often  called  the 
entropy  source  strength  -  as  a  unit  of  irreversibility: 


a 


d.S 

1 

dt 


(1-3) 


Combined  with  (1-2),  the  equation  expressing  the  first  law 
of  thermodynamics. 


dU  =  dQ  +  dW  ,  (1-4) 

becomes 

dW=dU  -TdS+T  d.S  =  d  W  +  T  d.S  (1-5) 

l  rev  l  v  ' 

where  dU  is  the  increase  in  internal  energy  of  the  system, 
and  (dW  -  d  W  )  Is  the  work  expended  or  the  energy 
dissipated.  The  general  form  of  the  rate  of  entropy  production  is 

a  =  2  X.  J.  (1-6) 

ii  '  ' 

i 

where  the  J\  T s  are  the  "fluxes"  denoting  the  speed  with 
which  the  irreversible  processes  occur,  and  the  ' s  are 
the  thermodynamic  "forces"  which  are  related  to  the  non¬ 
uniformity  of  the  system  (the  gradient  of  the  temperature  for 
example ) . 
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Since,  at  equilibrium  both  the  fluxes  and  the  forces 
disappear,  the  fluxes  are  surely  functions  of  the  forces.  For 
small  departures  from  equilibrium,  as  a  first  approximation 


Ji  "  X  LikXkJ 

where 

dJ, 

T  —  _ 

lk  “  3X^ 


(1-7) 


(1-8) 


Therefore 


a 


Z 

i,k 


L. ,  X.  X, 
lk  i  k 


(1-9) 


The  coefficients  are  called  the  phenomenological 

coefficients.  Since  (1-9)  is  always  positive  by  virtue  of 
the  second  law  of  thermodynamics,  the  matrix  [L^]  is 
positive  definite  and  is  symmetric,  i.e.,  . 

The  standard  procedure  to  find  the  phenomenological 
coefficients,  i.e.,  to  set  up  equation  (1-8),  is  summarised 
as  follows : 

From  the  rate  of  change  of  the  entropy,  which  is  determined 
by  the  variations  in  the  properties  of  the  system,  the 
entropy  source  strength  o  can  be  calculated.  From  (1-9)  if 
is  readily  noted  that  a  is  a  sum  of  scalar  products  of  the  form 


8  J. 

X.X,  . 
i  k 


- 


, 
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In  order  to  determine  the  rate  of  change  of  the 
entropy  of  a  system,  the  state  of  the  system  must  be  defined 
in  terms  of  certain  state  variables.  These  state  variables 
are  of  two  kinds:  l)  external  (e.g.  pressure,  temperature) 
variables  which  need  no  further  comment  and,  2)  "chemical" 
variables  if  reactions  occur  which  change  the  material 
composition  of  the  system.  Here  the  term  "chemical"  is 
being  used  rather  loosely;  excitations  of  internal  degrees  of 
freedom  constitute,  in  this  sense,  chemical  reactions. 

Thus  the  chemical  variables  need  to  be  examined  more  carefully. 
When  chemical  reactions  occur  in  a  system,  they  can  be 
described  in  terms  of  two  sets  of  thermodynamically  conjugate 
variables,  called  the  progress  variables  and  the  affinities. 

The  progress  variables  are  best  understood  as  follows 
Suppose  there  are  r  possible  reactions  in  a  system  containing 
N  substances.  The  reaction  equations  are  given  by 


/v 


c=/ 


£  L^lj]  (i-io) 


where  NL  is  the  mass  of  substance  i,  \JLr^  and  are  the 

left  and  right  stoichiometric  coefficients  of  substance  i  in 

the  reaction  a.  When  the  reaction  a  proceeds  by  a  differ- 

a 

ential  amount,  d  n^  the  change  in  the  number  of  moles  of  sub¬ 


stance  i  is  given  by 


cl/vu, 


I  c7(- 

ol/vu^ 


=.  d  a 


(1-11) 


Ay, 


IcK 


A  vi 


where  Av*  =  V/*  -  .  The  total  change  d 

becomes 


A- 

cIm  <■  =  Z1  A  v^o/£,^ 

<A~! 


(1-12) 


d  £>q  represents  the  differential  amount  by  which  reaction  a 
has  advanced.  The  parameter  is  known  as  the  progress 

variable.  When  the  reaction  has  proceeded  to  completion  to 
the  right,  =  1*  i.e.  A  molecules  of  substance  i 

have  been  produced. 

As  for  the  affinities,  well  before  a  practical 
measure  of  these  somewhat  "anthropomorphic  quantities”  was 
found  (Joos),  the  concept  of  chemical  affinity  had  been  widely 
used.  Van’t  Hoff  first  recognised  that  for  reactions  at 
constant  V  and  T,  the  free  energy  F  is  the  proper  measure 
of  affinity,  while  for  reactions  at  constant  p  and  T, 
the  Gibb's  potential  G  serves  this  purpose.  These  quantitie 
are  minimum  at  equilibrium;  the  greater  the  values  of  these 
functions  compared  with  the  minimum  values,  the  further  the 
system  is  from  equilibrium,  and  consequently,  the  more 
readily  the  reaction  will  occur. 


' 


>: 
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If  only  expansion  work  is  performed  on  a  system 


all  reactions  then  being  fixed,  the  process  is  reversible. 

The  quantity  -pdV  is  then  the  contribution  of  the  external 

variables  to  dW  .  If  however,  reactions  occur  at  constant 

rev 

volume,  the  process  is  not  reversible.  Therefore  other 
terms  must  be  subtracted  from  -pdV  in  order  to  obtain 
dW  .  These  terms  will  be  of  the  form 


rev 


2  Aa  d£, 
a 


where  the  A^’s  are  the  affinities,  and  the^^’s  are  the 


progress  variables.  Therefore 


The  first  law  of  thermodynamics  thus  takes  the  form 


=  T  dS-pdV  -  2  Aq  d£a,  (1-13) 


a 


from  which  we  obtain  an  expression  for  Aa: 


(1-14) 
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Just  as  -p  and  T  are  thermodynamically  conjugate  to 
V  and  S  respectively,  -A  is  thermodynamically  conjugate 
to  £4. 

Since  the  various  thermodynamic  potentials  are 
defined  by  the  equations 

G  =  U  -  TS  +  pV 
F  =  U  -  TS 

and  H  =  U  +  pV, 

the  following  expressions  for  A  also  hold: 


where  G,  F,  H  are  respectively  the  Gibbs  free  enthalpy,  the 
Helmholtz  free  energy  and  the  enthalpy. 

The  physical  variables  x  and  y  will  now  be 
introduced.  The  quantity  y  stands  for  either  S  or  T, 
while  x  stands  for  either  V  or  -p.  The  conjugate  variables 
will  be  expressed  by  X  and  Y.  Finally  P  denotes  the 
thermodynamic  potential  corresponding  to  x  and  y. 
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Equation  (1-13)  takes  the  form:  (x=V,  X=p, 

y=s,  Y=T) 


dP  =  -  X  dx  I  Y  dy  -  Z  d  ^ 


d-16) 


while  equation  ( 1 —15)  becomes 


A 


=  -(22 
\d£.. 


(1-17) 


Differentiating  with  respect  to  x  yields: 


3Ah|  =  /  a8  P 


12L 


(1-18) 


A  few  other  useful  formulae  can  be  derived,  and  are  given 
below. 


3  /A  ^ 
3  /j> 


!j 


3  A<A  /p) 


(l-19a) 


3  £, 


q_AA  /  3 


ax 


(l-19b) 


*  /Xi|  \"31t  l#<j  V3X  4y  '^r  AV’j 

The  fundamentals  of  the  thermodynamics  of  irrever¬ 


sible  processes  have  been  summarized  in  this  section, 
following  the  treatment  by  Bauer  (1965)* 
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2.  Classical  Sound  Absorption 

Using  the  principles  stated  in  the  preceding  section, 
the  rate  of  entropy  production  will  be  determined,  showing 
that,  in  the  case  of  an  isotropic  fluid  without  reactions. 


while  their  respective  conjugated  forces  are  the  gradient  of 
the  reciprocal  temperature,  and  the  flow  velocities.  The 
phenomenological  equations  will  then  be  set  up,  which  together 
with  the  equations  of  continuity  for  mass,  momentum  and 
energy,  constitute  a  system  of  homogeneous  equations.  When 
this  system  of  equations  is  linearized  and  solved,  the  complex 
wave  vector  k  for  the  sound  propagation  can  be  found.  A 
detailed  exposition  for  finding  the  wave  vector  k  now  follows. 


By  the  principle  of  conservation  of  energy,  the 


total  energy  content  within  an  arbitrary  volume  V  in  the 
system  can  only  change  if  energy  flows  into  or  out  of  the 
volume  through  its  surface.  With  the  help  of  Gauss T  theorem, 
the  differential  form  of  the  law  of  conservation  of  energy  is 
obtained : 


(1-25) 


at 


where  P  is  the  density,  e  is  the  total  energy  per  unit  mass 


and  the  total  energy  flux  per  unit  surface  and  unit  time. 

Here  e  includes  all  forms  of  energy  in  the  system  and  may  be 

written  in  the  form. 
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e  =  75-  v2  +  Uj  (1-26) 

1  2 

where  75-  v  is  the  specific  kinetic  energy,  and  u  is  the 
specific  internal  energy.  Similarly  includes  all  forms 

of  energy  fluxes:  a  convective  term  p  ev,  an  energy  flux 
Tl*v  due  to  the  mechanical  work  performed  on  the  system, 
and  finally,  the  heat  flux  T  : 

Je  =  pev  +TT-V  +  .  (1-27) 

The  conservation  law  for  energy  of  a  fluid  in  the  absence 
of  external  forces  can  thus  be  written: 

pf|=  -  V-  m-v  +  ?q)  (1-28) 

where 

it  -  It  +  7  (1-29) 

is  the  so-called  differential  following  the  motion. 

Remembering  that  in  a  fluid,  only  differences  of 
pressure  give  rise  to  a  force  (constant,  normal  or  shearing 
stresses  do  not  produce  a  force)  the  conservation  law  of 


momentum  can  be  written  as 
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P  fir =  -  ^  •  (1-30) 

Combining  equations  (1-26)  and  (l-28)  gives  the  balance  of 
internal  energy. 

PH¥  =  7  fj  d-31) 

It  is  now  useful  to  split  the  pressure  tensor 
into  two  parts:  the  hydrostatic  pressure  p,  and  the  symmetric 
viscous  pressure  tensor  TT  This  is  justified,  since  the 
components  of  the  stress  tensor  are  the  following: 


P 


-  a  it. 


( 1 -32a) 


n. 


( 


( 1 -32b) 


where 


x,y,z 


v. 

l 


and  \\  is  the  shear  viscosity.  For  a  monatomic  gas,  the 
bulk  viscosity  is  identically  equal  to  zero.  The  hydrostatic 


pressure 


St-  <x 
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p  =  j  (rrxx  +vyy  +ttzz). 


(1-33) 


The  balance  equation  for  the  internal  energy  and  the  law  of 
conservation  of  mass  become  respectively: 


P  =  -  P  V  •  V*  -  Us  :  V  v>  -  V  •  ,  (1-34) 


and 


q 


(Vsp-1)  (1-35) 


From  the  last  two  equations. 


Kir  +  >,$f)s'TT*:V7?~  v'Jt  (1'36) 


which  can  be  inserted  into  the  Gibbs  relation 


T 


dS  =  du  u  d\/ 
dt  dt  \  dt 


(1-37) 


Thus  the  entropy  balance  equation  is  obtained: 


pf=-V-T  +  <r 


(1-38) 


with  the  entropy  flow 


J 


J 


q 


s  T 


(1-39) 


f 

■•  V  ,tt-«  ■  . 

,  ;  -r  -C  -V  3c:"*  .. 

9  -  ::  . 
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and  the  entropy  production 


<r  =  -  X  •  V (T)  ~  TP  :Vv» 

T*  T 


>  o 


(1-40) 


As  mentioned  before,  the  entropy  production  is  a  sum  of 

products  of  fluxes  and  thermodynamic  forces,  the  former 

being  related  to  the  latter  through  the  phenomenological 

equations  (1-7).  Since  the  phenomenological  coefficients 

can  be  temperature  dependent,  new  coefficients  can  now  be 

introduced  into  which  such  factors  as  ^  and  can  be 

T 

absorbed . 


The  familiar  expression 


Jq  =  -  P  V  T> 


(1-41) 


where  (p  is  the  heat  conductivity,  is  obtained  through  such 
simple  considerations,  from  equation  (1-40).  Similarly, 
equations  (1-32)  yield  in  turn 


TTS  =  -2n(vX  -  J_\7-X  U) 


(1-42) 


where  U  is  a  unit  tensor. 


I 
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The  present  treatment  follows  in  broad  outline  and 

p 

is  a  simplification  of  the  treatment  by  De  Groot  and  Mazur  , 
but  will  henceforth  be  restricted  to  one  dimension.  In 
such  a  case,,  the  laws  of  conservation  of  mass  and  momentum  are 

do...  ,  p  d.y» 

C\  T  1  c)  /£ 

and 

p  dvx__  c)  E 

>  d t  rr 


=  0  (1-^3) 

4-4  4-4  =  O  (1-44) 

■i  c)  # 


Again,  restricting  the  problem  to  one -dimension,  and  inserting 
the  phenomenological  equations  (l-4l),  (l-42)  into  the 
energy  law  (l-3?),  the  following  expression  is  obtained: 


pr +p 


(1-45) 


Now  it  is  well  known  that  for  a  monatomic  gas  whose 
molecules  possess  translational  kinetic  energy  only,  the 
internal  energy  of  the  gas  equals  75-  kT.  Moreover,  from  the 
equation  of  state  of  the  gas  (k  is  Boltzmann's  constant) 


(1-46) 


>9  !•  l^oi-nsaoi. 


1 
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where  m  is  the  mass  of  a  molecule  of  the  gas.  Therefore 
the  internal  energy  per  unit  mass  becomes 

U  =  |  £  •  (1-47) 

The  three  preceding  equations  have  to  be  linearized,, 
in  order  to  handle  them  easily  mathematically.  In  the  case 
of  moderately  weak  ultrasonic  waves ,  the  usual  procedure  for 
linearizing  these  equations  is  to  express  p  and  p  as 

o  =po  (  1  +  £  )  ,  p  .  po  (  1  +  <5  )  .  (1-48) 

Here  ^  q  and  pQ  are  respectively  the  density  and  the 
pressure  of  the  fluid  in  the  absence  of  perturbation.  £  and  £ 
are  small  compared  to  unity,  while  V*  the  flow  velocity 
is  very  much  smaller  than  the  sound  velocity.  Therefore, 
powers  higher  than  one  and  products  of  these  quantities  or  of 
their  derivatives,  can  be  neglected.  Substituting  for  ,  p 
and  u,  equations  (1-43) *  (1-44)  and  (1-45)  become: 


de  ,  c)  v* 
3  t  3  ^ 


=  0 


(1-49) 


8  Vy.  ,  JaL  -4-4-  -  o  and  (1-50) 

Q  £  i°0  drp  1  P° 


4  +  2^-  -  -Sr-  ^ 4(<5_£)  -  0  ■  (1-51) 

3  or  3  p„  k  3^ 


edr  nsritf  *isJ 


H  *  j^LC: 

vVo 
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Eliminating  v  from  the  above  equations ,  there  results 

-X. 


.J-Lr  C)1  £ 

p0  d ^ 

and 


±  A_ 


=  o 


(1-52) 


a J  f  be 

b  t  3  ~5T 


(1-53) 


Equations  (1-52)  and  (l-53)  are  the  equations 
governing  the  propagation  of  sound  waves  of  small  amplitude 
in  a  gas  treated  as  a  continuous  medium.  Assuming  that 
the  sound  wave  is  propagated  as  a  damped  progressive  wave, 
e  and  5  must  be  of  the  form 

i  ( \{x  -  u>t )  ,  R  R  i  ( Kx  -  ^t ) 

e  =  e  e  v  x  '  and  6  =  6  e  VIX  '  . 

o  o 


Substituting  these  values  into  equations  (l-52),  (l-53)  yields 
two  homogeneous  equations  for  and  5  . 

( uj  +  l  i_  Jx.  uJ  K  \  -  O  (1-54) 

V  i  J  p, 

€'  /x  +t  it  EL  \  +  Sc  (  -  L  JL<Ll£  (1-55) 

l  3  3  kuu  J  \  3(°0k  J 


“ 


O  r 


19 


Since  the  determinant  of  these  equations  must  be  zero  for 


consistency,  we  obtain  a  quadratic  equation  in  K  ,  namely 


K 


-  L  Jl  w. 

+  (  2. 

+  N' 

.  i  p.n 

3 

fa 

X  fa.  +  L  UL-ff 
”  fa  3  fa  ^ 


+  i  JL  Mr- 

3  p0  J 


+  UJ 


*  O 


(1-56) 


An  approximation  to  the  classical  sound  absorption 

coefficient  o{  ,  ,  which  is  often  encountered  in  the 

' class 

literature  can  be  derived  from  equation  (1-56)  by  means  of 
various  substitutions.  Firstly 


V 


r  P  1/2 
(1  12.  ) 

9  Po  } 


is  the  adiabatic  sound  velocity  at  low  frequencies  for  a 

K  K  V- 

monatomic  gas.  Next  let 


k4  =  p~  = 

1  Po 


o 


CD 


<f  =  15  k 


Further,  following  Greenspan 


and 

3 


^  2, 

let  R  =  -fa  Vq  where  R  is  the  Reynold's  number  taken  from 

t]  UJ 

the  field  of  gas  dynamics.  Equation  (1-56)  becomes 


1  i_  -  lo  1 
U  3  R* 


+  fa-  I  +  Xii 


6  K  J 


+ 


=  O 


or 


ki  =  ft 


^  xa) 


'A 


^  -UK 


Thus  fa  =  (.  +-ZJ_  -  J !_£i  fa 
1  4  7\  72 


- 

V3  a 


/ 

R3 


-f 


(1-57) 


which  is  obtained  from  the  preceding  equation  by  series 


approximation . 


f  xo'K 
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In  the  zeroth  approximation,  the  velocity  is 
constant,  and  there  is  no  attenuation.  The  first  approxima¬ 
tion  gives  the  "classical"  attenuation  per  unit  length 


<A  class 


P 


o 


class  A  o 

2  TT 


7  1 

"Fr 


(1-58) 


where  is  the  wave-length  at  low  frequencies.  Equation 

(1-58)  can  be  shown  to  be  equivalent  to  the  usual  approximation 


(1-59) 


Here  is  the  specific  heat  at  constant  pressure,  y  is 

the  ratio  of  the  specific  heats  and  f  is  the  number  of 
vibrations  per  unit  time  of  the  sound  wave. 

Indeed  the  equivalence  of  (1-58),  (1-59)  follows 
5 

easily,  since  y  =  in  the  case  of  a  monatomic  gas  and 

A  =  7  by  the  kinetic  theory  of  gases,  where  Qv  is  the 

specific  heat  at  constant  volume.  Thus  the  classical  sound 
absorption  is  seen  from  (1-59)*  to  be  directly  proportional 
to  the  square  of  the  frequency,  and  inversely  proportional 


to  the  pressure. 


4 
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3.  Relaxational  Absorption 

When  a  fluid  is  given  a  periodic  deformation,  as 
is  the  case  when  a  sound  wave  is  transmitted  through  the 
fluid,  a  phase  lag  exists  between  the  stresses  and  the 
strains,  which  is  responsible  for  the  absorption  of  sound. 

In  monatomic  gases,  this  phase  lag  is  due  entirely  to  vis¬ 
cosity  and  heat  conduction.  In  diatomic  and  polyatomic  gases 
however,  the  absorption  of  sound  is  also  due  to  the  failure 
of  the  internal  degrees  of  freedom  to  remain  in  thermal 
equilibrium  with  the  external  degrees  of  freedom.  The 
process  by  which  the  temperature  of  the  internal  degrees  of 
freedom  adjusts  itself  to  the  temperature  of  the  external 
degrees  of  freedom  is  a  typical  relaxation  process.  Due 
to  the  passage  of  the  sound  wave,  the  "translational  tempera¬ 
ture”  is  increased.  In  order  to  be  again  in  equilibrium  with 
the  new  translational  energy,  the  internal  energy  must  also 
increase;  i.e.  energy  must  flow  from  the  translational  to 
the  internal  degrees  of  freedom.  This  can  only  occur  through 
collisions,  and  takes  a  finite  time.  The  time  needed  to 
decrease  the  maladjustment  by  a  factor  e  ^  is  called  the 
relaxation  time.  Such  a  process  is  called  a  reaction  in  the 
language  of  the  thermodynamics  of  irreversible  processes  and 
the  value  of  the  relaxation  time  depends  on  whether  the  reaction 
proceeds  isothermally  or  adiabati cally . 


ctt  bos  a»aeei38  aeswctdd  arfalxs  gf  C  ©a  iriq  «r  J  ^ 
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Again,  following  Bauer  (1965)*  the  phenomenological 
equations  are  easily  set  up  as  follows.  Combining  equations 
(1-3) *  (l-5)  and  (1-13)  yields  the  rate  of  entropy  production. 

After  absorbing  the  factor  ^  into  the  phenomeno¬ 
logical  coefficients,  the  phenomenological  equations  become: 


2 

(3 


L 


A, 


(<*=  1,  .  .  .*  r)  .  (1-60) 


In  the  case  of  a  single  relaxation  process,  equation 
(l-6o)  takes  the  form: 


♦ 

^  =  LA  .  (l-6l) 

L  may  be  assumed  to  be  constant  within  a  small  region. 
However  A  depends  on  the  variables  x(  =  V  or  -p), 
y(:  S  or  T)  and  £>  .  Within  a  small  region  about  the 
point  (xo,  y0,£,o)-- 

A  =A(^y„y0')  (<£-s,) 

v  ^ 


(1-62) 
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If  the  point  (x  ,  y  , £  )  is  an  equilibrium 

o  o^o 

point,  A(x  ,  yQ)£0)  vanishes*  since  as  stated  earlier,  the 
affinities  vanish  at  equilibrium.  From  equations  (l-6o), 
(l-6l)  the  so-called  relaxation  equation  (1-62)  is  obtained. 


Moreover  by  holding  the  external  variables  x  and  y  constant 
at  x  =  x  ,  y  =  yQ  the  last  equation  becomes  a  linear 
homogeneous  differential  equation  whose  solution  is 


£o)  “  (^>  t-.o  ]  (1-64) 

with  the  reciprocal  relaxation  time 

\  =  "  L(|f)  *  (1-65) 

L  xy  xy 

If  the  experiment  is  performed  under  different 
conditions,  so  that  X  and  y  are  kept  constant,  rather 
than  x  and  y,  a  different  relaxation  time  will 

result.  The  relation  between  these  various  relaxation  times 
is  given  by  equations  (l-65)  and  (l-19b) 


* 
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t 


t 


*3 


=  (M  (Ml 
^  ©5T/,3 


represents  the  differentiation 


the  operator  / _ ) 

l  dp  Ja -j 

carried  out  during  equilibrium,  while  the  operator 


^4 


represents  the  differentiation  carried  out  during  an  arrested 
reaction.  The  values  obtained  by  means  of  these  operators 
are  often  referred  to  as  quasi-static  and  instantaneous  or 
frozen-in  values  respectively. 

For  the  four  variables  p,  V,  T  and  S  ,  the 
various  relaxation  times  are  related  as  follows: 


(1-66) 


;  (1-67) 

;  (1-68) 
•  (1-69) 


K  and  K  are  the  isothermal  and  the  adiabatic  compress!- 
I 

bilities  respectively.  The  indices  o  and  °o  refer  to  the 
quasi -static  and  instantaneous  values  respectively. 


rs9fl  ■ 


baa  Isnnsrfjo: 
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Unfortunately,  the  progress  variables  cannot  be 
observed  themselves,  and  the  relaxation  times  cannot  be 
determined  as  easily  as  equation  (1-64)  seems  to  indicate. 
They  can  only  be  found  from  the  behaviour  of  another  depen¬ 
dent  variable  z,  which  is  a  function  of  x,  y  and  <E^  . 
Therefore,  if  instead  of  being  kept  constant,  x  and  y 
are  allowed  to  vary,  with  the  help  of  (1-65)  and  (l-19a), 
equation  (l-63)  becomes: 


On  the  other  hand  z  =  z (x,y,£j  can  be  expanded  linearly 
about  x  ,y  as  follows: 


z,z  can  be  eliminated  from  the  last  two  equations,  yielding 


2  6 


From  equation  (l-19b),  it  follows  that 


and  that 


a/y 


m 

'  J  'a# 


=pi) 


J2± 

\3* 


<3 


Equation  (1-72)  now  becomes: 


l  +  L  spy  jL- 


d  t 


(ri°)  =fe+(§£)  *«•»  &■](*-*) 


+ 


alj  t  >pu  j2l_ 

^3  if  dt. 


(1-73) 

(ft)  • 


Assuming  that  the  state  variables  are  of  the  form: 


(x  -  x  )  =  Cef^+f 
v  o' 


Q 

the  differential  operator  is  simply j  tU  .  Thus  equation 

(1-73)  takes  on  the  form 
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or 


(*“>)  =(^  +(^)J 


'K 


; 


where 


3  3 

3  rf 


and 


ii 


ell 


4 


1U 


=  ( hi 

3  Ajj 


i  +  j  w  r 


a>- 


tj  /  -+  j  dJ  L  ^ 


A 


1  4"  J 

/  -f  j  uj  t#  j 


The  derivatives 


] 

(3  Wy 


and 


ajY^ 

^  l 


(l-74a) 


(l-74b) 


are  clearly 


the  effective  state  coefficients,  which  take  into  account 
the  finite  reaction  rates.  These  equations  can  now  be 
applied  to  the  case  of  acoustic  relaxation. 


If  we  neglect  the  transport  phenomena,  the  adiabatic 

character  of  the  process  is  guaranteed.  In  the  absence  of 

the  dissipative  factors,  the  equations  of  continuity  for 

momentum,  mass  and  energy  yield  the  Laplacian  formula  for 

the  wave-vector  k  : 

o 


k 

o 


2 


(1-75) 
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where  Kg°  is  the  adiabatic  compressibility,  in  the  absence 

of  relaxation.  If  the  internal  transitions  are  now  taken 

into  account,  the  adiabatic  compressibility  must  be  replaced 

eff  eff 

by  an  effective  adiabatic  compressibility  can 

be  easily  found  in  terms  of  £  and  £  .  To  do  this  the 

°  yz  u  xy 

following  choice  of  variables  is  made: 


x  5  -  p. 


Using  equation  (1-74)  with 


^  5  r 


a  "normalized"  effective  adiabatic  compressibility  is  obtained 


K, 


eff 


K, 


o 


.  /j< 

I  -f  J  uj  L  vs. 

I  +j  UJ  5 


The  locus  of  KQeff/K  0  can  be  obtained  by  means  of  the 

conformal  mapping  W  =  /)  4  j  ^ vA.  onto  the  positive  real  axis 

.  )  +  i  I  £  s 

z  =  l/o).  This  gives  a  semi -circle  of  diameter  €  =  /  _ — 

€  is  called  the  relaxation  strength.  The  maximum  of  the 


imaginary  part  occurs  when  UJ  =  - 


r 


This  locus. 


which  is  known  as  the  Cole  plot,  is  useful  in  determining 


a  I  '  > 


■ 


■  -  .  .  ■  ' 
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whether  sound  absorption  or  dispersion  is  due  to  a  single 
relaxation  process  or  to  a  multiple  relaxation  process:  the 
curve  is  a  semi-circle  in  the  former  case  only. 

Under  these  conditions,  equation  (1-75)  must  be 
replaced  by: 


k 


=  K, 


eff 


CD 


(1-76) 


where  kr  is  the  complex  wave  vector 


k 


(JD 


r  ~  C 


-  J 


* 


(1-77) 


If  the  value  of  kr  given  by  (1-77)  is  substituted  into 


(I-76),  and  both  sides  are  multiplied  by  C 
an  obvious  simplification  yields 


2 


I 


o 


pK 


C 

c1 


_  _  i  a c?  =  K 

■a  J  cue  W 


U J 


1  ,  A  *\ 

I  ~t  LL)  i 


I  +  t\ 


_  j  u>  (i  r,S _ LinX 


P* 


I  +  uT  L 


(1-78) 


The  real  and  imaginary  parts  of  (1-78)  will  be  equated 
separately  in  what  follows.  Comparing  first  the  real  part. 


c 


£L_  — 


c 


)  _  3. 

uj3. 


1  +  u  tr 

1  +  u)1  tq 


(1-79) 
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>  * 

Now,,  CQ  is  not  very  different  from  C.  Thus  d  £- „ 


—  t A 


A 


M  \ 

,  where  A 


also  be  written  as 

cu*  V7T 

wavelength  of  the  sound  wave.  Moreover,  in  most  cases 

/  j  v  A  / 

id  X)  /  it1  <3<  /,  Under  these  conditions,  (l~79)  can  be 
approximated  by 


C  2  i  _|_  a)2  f*  T 

o  CVS  lpS 

- - - - 


p  - 

l  +  or  -f- 


pS 


However,  from  (1-66) 


oo 


f  t  =  t2  —  =  (i  _  € ) 

LVSLpS  pS  K  o  1  pS  u  >  ’ 


so  that  finally 


C 


2 


C 


2  2 

£  or  f 


=  1  - 


pS 


1  + 


pS 


Plotted  on  a  logarithmic  frequency  scale,  this  function 
exhibits  a  dispersion  step  with  an  inflection  point  at 


0) 


pS 


If  on  the  other  hand,  the  imaginary  parts  of 
(1-78)  are  compared,  it  is  easily  seen  that 

A  1 u  (Us  -  4) 

is 


UJ  C 


 * 


A  Q. 

Jr  C 


/+  cu*  ?*< 


may 
is  the 
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o<  A  is  called  the  molecular  absorption  index,  and  is 
designated  by  ju  or  .  Using  (l-8o)  a  little  simplification 
yields : 


IT  G 


(1-81) 


and 


u  .  cA  )\  -  TT  € 

7 


uJ 


a, 

/+  4  4 


(1-82) 


I  A 

When  plotted  against  In  co  ,  both  Cc  and  yx  are  bell- 
shaped  curves.  The  maxima  of  these  curves  are  at  = 


and 


u j  = 


\r 


i 


f- 


respectively. 


PS 


C 


The  curves  of 


C 


//  c 

Li  o 

’  — U- 

C 


and 


appear  m 


fig.  (l)  a),  b)  and  c)  respectively. 


Acoustic  measurements  can  therefore  yield  the 

values  of  ^  and  *  Equation  (l-8l)  in  particular 

constitutes  the  basis  for  the  measurements  reported  in  this 

thesis;  the  molecular  absorption  index  is  thus  seen  to  be 

related  to  the  two  adiabatic  relaxation  time  T.  0  and  T1Ta . 

pS  VS 

On  the  other  hand  the  relaxation  times  can  be 
calculated  by  means  of  equations  (1-65) •  However,  the  thermo¬ 
dynamics  of  irreversible  processes  does  not  provide  any  information 
about  the  progress  variables  or  about  the  affinities.  Reaction 


' 
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Fig.  1 


a)  (CQ/C)2  vs  In  co 

b)  [x  (CQ/C)2  vs  In 

c) 


l_i  vs  In 


In 
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kinetics  must  therefore  be  applied  in  order  to  evaluate 


these  parameters.  For  the  sake  of  simplicity  and  following 
4 

Kneser  C  will  be  calculated  only  for  the  special  case 


where  the  reaction  equation  has  the  following  form: 


A  -f- 


A  +-  A 

n - —  o 


(1-83) 


This  equation  describes  the  de-excitation  and 

excitation  of  internal  degrees  of  freedom  by  binary  collisions. 

A  denotes  the  unexcited  molecule  and  A  represents  the 
o 

collision  partner  (A  or  A  )  which  remains  unchanged 
during  the  collision.  Although  equation  (1-83)  is  only 
one  of  the  several  reactions  which  are  possible  in  principle,, 
it  is  the  only  reaction  which  occurs  as  the  result  of  colli¬ 
sions  in  a  two  level  gas. 

The  population  of  the  ground  state  and  that  of  the 
excited  state  change.,  during  the  reaction  to  the  right, 
according  to  the  equation: 


d  n>  d  n* 

_ o  _  _  _ 1 

dt  dt 


Rn 


1 


where  nQ  and  n^  are  the  number  of  molelecules  per  unit- 
volume  in  the  ground  and  in  the  excited  states  respectively; 
R  is  the  rate  constant  to  the  right. 

A  similar  equation  describes  the  reaction  to 


the  left : 
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R'n 

o 


Here  R’  is  the  rate  constant  to  the  left.  If  the  molar 


number  n  is  introduced  as  the  progress  variable  £, 
0  .  d  "n^  d  1n“ 

the  time  derivative  includes  —  -  and 


dt 


dt 


can  be  expressed  by 


4  =  Rn-L 


R'n  —  Rn, 
o  1 


1 


(1-84) 


At  equilibrium.,  =  0  and  n0o^nio  =  where  the 

additional  subscript  o  (on  the  n's)  denotes  the  equilibrium 
values.  At  small  displacements 


n 


1 


n 


oo 


n 


lo 


+  A  no 

+  A  n! 


n  -+-(£-£)  n 

OQ  v<->  _  oo 

nT  1  (>  -  jr  T"  ~  n 
lo  ^  y  lo 


1  +■ 


oo 


This  value  substituted  into  (1-84)  yields 

i  =  -  =  -  r  >)  k-s0)  •  u-85) 

00  lo  oo 


When  the  experiment  is  conducted  under  isothermal  conditions,, 

?  remains  constant  during  the  relaxation  process .  The 
*“■>0 

solution  of  (1-85)  is  then 


.  ~~jr  .  o  ~  3‘.  O'  :  ijti  £  r\  •  *.J  13b  li  ’  o  ;  > 
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t=  o 


e 


with  an  isothermal  relaxation  time  given  by 


zy1  =  R  ( 1  +-  )  =  (H  t-H1  )  .  (1-86) 

oo 

When  the  reaction  represents  the  excitation  or  the  de¬ 
excitation  of  a  molecule.,  the  rate  constants  R  and  R’  are 
simply  proportional  to  the  probabilities  f  and  f^Q: 


R  = 


R'.-r1 

c 


where  is  the  mean  collision  time;  f^Q  is  the  probability 
that  an  excited  molecule  will  return  to  the  ground  state  when 
colliding  with  another  one,  while  fQ^  is  the  probability  that 
a  molecule  in  the  ground  state  will  be  excited  when  colliding 
with  another  molecule.  Equation  (l-86)  shows  that  the  iso¬ 
thermal  relaxation  time  is  inversely  proportional  to  the  sum 
of  the  rate  constants. 

The  experimental  values  of  the  sound  absorption  or 

of  the  velocity  dispersion  yield  the  adiabatic  relaxation 

times  T  ~  and  urr1  ,  while  reaction  kinetics  calculations 
pS  VS 

or  •  I  ft 

order  to  compare  the  experimental  values  with  the  theoretical 


yield  the  isothermal  relaxation  times 


t 

pT 


o 
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values,,  the  former  must  be  converted  into  Isothermal  values 
by  means  of  equations  (I-67)  and  ( 1 —68 ) .  The  difference 
between  the  isothermal  and  the  adiabatic  relaxation  times  is 
often  minimal,  especially  in  the  case  of  vibrational  relaxation. 
It  is  not  so  however  in  the  case  of  rotational  relaxation. 

The  various  relaxation  times  are  clearly  indicated  on  fig.  (2). 

The  following  inequalities  always  hold: 
t  pt  >  fpS  ^  /VS 

This  distinction  between  the  various  relaxation 
times  has  often  been  ignored  in  the  literature.  Nowhere  has 

it  been  expressed  more  clearly  however,  than  in  the  papers 

1  4 

by  Bauer  and  by  Kneser. 


4.  Relaxational  Absorption  in  Hydrogen. 

In  order  to  compare  experimental  results  with 
theoretical  calculations,  it  is  best  to  effect  measurements 
on  a  simple  molecule.  Limiting  ourselves  to  the  problem  of 
rotational  energy  transfer,  the  hydrogen  molecule  at  low 
temperature  constitutes  such  a  simple  model  for  the  following 


reasons : 


' 
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Fig .  2 


a) 

b) 

c) 
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Tg :  external  temperature 

1\:  internal  temperature 
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l)  depending  on  the  isotopic  constitution  of  the  molecules,, 
only  two  levels  are  important  at  temperatures  ranging  from 
20  to  150°  K  --  when  more  than  two  levels  are  occupied, 
multiple  relaxation  processes  occur,  which  are  rather  diffi¬ 
cult  to  analyse;  2)  due  to  the  very  small  moment  of  inertia 
of  the  hydrogen  molecule,  the  rotational  energy  levels  are 
widely  spaced;  3)  the  relaxation  strength  e  can  be 
computed  easily  since  physical  properties  such  as  the  specific 
heats  are  well  known  in  the  case  of  hydrogen.  These  three 
points  will  now  be  examined  more  closely. 

The  hydrogen  molecule  is  essentially  a  rigid 
rotator  whose  energy  levels  are  given  by  the  usual  formula 

~2 

Ej  =  J(J  +  l)|j- 

where  I  is  the  moment  of  inertia  of  the  rotator,  J  is  the 
rotational  quantum  number,  ii  is  Planck Ts  constant  divided 
by  2t r  .  Depending  on  the  orientation  of  the  nuclear  spins 
of  the  two  hydrogen  atoms  which  constitute  the  molecule, 
hydrogen  can  exist  in  two  modifications:  orthohydrogen  when 

the  total  nuclear  spin  is  one,  and  parahydrogen  when  the  total 
nuclear  spin  is  zero.  Quantum  mechanical  considerations  show 
that  only  the  odd-numbered  rotational  states  can  be  occupied  in 
the  case  of  orthohydrogen,  and  only  even-numbered  rotational  levels 
can  be  occupied  in  the  case  of  parahydrogen. 
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The  thermodynamical  equilibrium  between  the  ortho 
and  the  para  modifications  is  given  by  Boltzmann's  distri¬ 
bution  law,  and  simple  calculations  show  that  at  room  tem¬ 
perature,  orthohydrogen  (o-H^)  represents  75^  of  the 
equilibrium  mixture.  At  20°K  however,  the  equilibrium  mixture 
is  over  99%  parahydrogen  (p-H2) .  Transitions  between  the 
two  modifications  do  not  occur  except  in  the  presence  of  a 
catalyst.  For  example,  in  the  absence  of  a  catalyst,  normal 
hydrogen  (75%>  ortho  -  25^  para)  does  not  convert  to  pure 
parahydrogen,  even  when  cooled  down  to  20°K.  The  relative 
concentration  of  ortho  and  para  molecules,  3:1  remains  the 
same.  As  a  result,  the  ortho  and  para  modifications  can 
be  treated  as  two  different  gases.  In  the  presence  of  a 
catalyst  however,  almost  pure  parahydrogen  can  be  obtained 
by  cooling  hydrogen  down  to  20°K. 

The  occupation  of  the  rotational  energy  levels 
as  a  function  of  temperature  for  p-H^  and  o-H^  appears  in 
figure  (3)  a)  and  b)  respectively.  As  can  be  seen  in  these 
figures,  below  100°K,  parahydrogen  is  a  two-level  gas,  which 
can  therefore  exhibit  a  single  relaxation  process,  under  the 
influence  of  a  sound  wave. 

The  fact  that  the  hydrogen  molecule  has  a  small 
moment  of  inertia  is  of  great  importance;  indeed,  the  spacing 
between  the  rotational  energy  levels  of  a  molecule  is 


4o 


Fig.  3  a) 


Energy  level  population  in  parahydrogen . 
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Fig.  3  b) 

Energy  level  population  in  orthohydrogen. 


n 


inversely  proportional  to  the  moment  of  inertia  of  the  mole¬ 
cule,,  and  transitions  are  less  likely  to  occur  when  the 
energy  levels  are  widely  separated.  Therefore,  one  can 
conclude  that  molecules  with  small  moments  of  inertia  have 
comparatively  long  relaxation  times.  In  fact,  most  molecules 
have  such  large  moments  of  inertia  that  their  rotational  relaxa¬ 
tion  times  are  too  short  to  be  measured  accurately:  hydrogen, 
deuterium  and  hydrogen  deuteride  are  a  few  exceptions. 


Equations  (l-8o),  (l-8l)  and  (1-82)  contain  the 
factor  e,  the  relaxation  strength.  As  mentioned  previously, 
e  is  the  diameter  of  the  "Cole  plot",  and  is  equal  to 
(l  -  the  help  of  equations  (1-66),  (1-67) 

and  (1-68) ,  and  putting  K^/K^  =  1,  one  can  write 


e  = 


0  oo  o©  n 

C1T  C  -  CTr  Cu 
V  p _ V  p 

0  Co 

V  CP 


Introducing  CT  as  the  rotational  specific  heat,  i.e.  that 


part  of  the  specific  heat  due  to  the  rotation  of  the  molecule. 


and 


Consequently, 


C 


oo 


p 


c 


0 


P 


-  c 


CO  0 

n  _  f  u  _  fi 


C T (c  0  -  c„°) 

v  p _ V  ' 

cv°  (CP°  ■  C,) 


The  relaxation  strength  is  thus  seen  to  be  directly  related  to 
the  values  of  the  specific  heats. 
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CHAPTER  II 


EXPERIMENT 


1.  Experimental  Methods 

The  methods  which  have  been  devised  for  measuring 
the  absorption  and  the  dispersion  of  sound  waves,,  can  be 
classified  into  three  groups:  l)  mechanical,  2)  optical,  and 
3)  electrical.  In  one  section  of  a  conprehensi ve  paper  on 
the  absorption  of  sound  in  fluids,  Markham,  Beyer  and  Lindsay 
made  a  critical  examination  of  these  methods,  and  the 
following  paragraphs  are  based  on  this  paper. 

Mechanical  methods  make  use  of  the  radiation 
pressure  which  is  the  difference,  due  to  the  sound  beam, 
between  the  pressure  on  one  side  of  a  rigid  wall,  and  the 
pressure  on  the  other  side  of  the  wall,  in  the  medium  at  rest. 
The  forces  being  measured  are  very  small,  and  often  of  the 
same  order  of  magnitude  as  other  extraneous  forces  such  as 
the  surface  tension.  Cavitation  effects  on  the  other  hand, 
set  an  upper  limit  to  the  sound  intensity.  Other  factors 
such  as  reflections  from  the  walls  of  the  container  and 
divergence  of  the  beam  at  lower  frequencies  limit  the  use¬ 
fulness  and  the  accuracy  of  this  method.  Above  10  Mc/sec,  the 
accuracy  is  ±  5/,  while  below  3  Mc/sec,  hardly  any  results 
obtained  by  this  method,  can  be  relied  upon. 
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Optical  methods  are  based  on  the  fact  that  when  a 
gas  is  subjected  to  a  sound  wave,  the  alternating  regions  of 
compression  and  rarefaction  act  as  diffracting  elements  for 
a  beam  of  light  at  right  angles  to  the  sound  path.  The 
sound  beam  transforms  the  gas  into  a  diffraction  grating. 

It  is  easy  to  see  that  such  a  method  imposes  stringent 
experimental  requirements:  alignment  for  instance  is  very 
critical.  Divergence  of  the  sound  beam  sets  a  lower  limit 
to  the  range  of  frequencies  at  which  experiments  can  be 
conducted.  Optical  methods,  like  mechanical  methods  are 
most  trustworthy  above  10  Mc/sec. 

The  majority  of  measurements  on  sound  absorption 
and  dispersion  in  gases  have  been  obtained  by  electrical 
methods,  which  include  all  methods  making  use  of  a  micro¬ 
phone.  Although  piezoelectric  transducers  are  the  commonest 
type  of  microphones,  transducers  of  the  condenser  type  with 
solid  dielectric  have  proved  to  be  equally  reliable  and  more 
versatile.  Electrical  methods  can  in  turn  be  classified  into 
three  groups:  a)  interferometric,  b)  direct,  and  c)  pulse. 

The  first  ultrasonic  interferometer  was  built  by 
Pierce  in  1925*  Shortly  after,  Hubbard  made  significant 
contributions  to  the  field  of  ultrasonic  interferometry. 

Most  of  the  earlier  results  on  the  absorption  and  dispersion 
of  sound  in  gases  were  obtained  by  this  method.  An  acoustic 
interferometer  is  essentially  an  electromechanical  system 
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composed  of  a  column  of  fluid  set  into  longitudinal  vibrations 
at  one  end  by  an  oscillating  plate,  and  terminated  at  the 
other  end  by  a  reflecting  plate.  Multiple  reflection  takes 
place  between  the  source  and  the  reflector;  as  the  length 
of  the  sound  path  is  varied,  cyclical  changes  of  phase  and 
amplitude  of  the  forced  vibrations  occur,  the  reflector 
passing  through  successive  positions  of  resonance  in  the 
fluid  medium.  Due  to  the  nature  of  the  transducer,  the 
variations  of  pressure  at  the  surface  of  the  vibrator  manifest 
themselves  as  variations  of  current,  or  sometimes  of  frequency, 
in  the  electric  circuit  coupled  to  the  transducer.  Since 
it  is  essentially  a  problem  of  standing  waves,  these  cyclical 
changes  occur  each  time  the  reflector  is  moved  by  half  a 
wavelength.  The  interferometer  thus  lends  itself  most 
conveniently  to  the  measurement  of  the  sound  velocity.  At 
first  glance,  it  would  appear  that  the  absorption  coefficient 
could  be  readily  obtained  from  the  ratio  of  the  amplitude  of 
two  successive  maxima.  In  practise  however,  many  factors 
such  as  the  reflection  coefficient  of  the  reflector  have  to 
be  taken  into  account  when  calculating  the  absorption 
coefficient . 

6 

Pumper  criticized  this  method  and  demonstrated 
that  the  interferometer  can  be  made  to  yield  differing  results 
by  slight  modifications  of  technique.  The  scatter  of  the 


. 
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results  obtained  by  this  method  for  the  values  of  the  sound 
absorption  coefficients  seem  to  justify  Pumper's  criticisms. 

As  stated  previously  however,  the  interferometer  remains 
most  useful  and  reliable  for  measuring  the  sound  dispersion. 

In  the  direct  method,  a  transducer  of  radius  a 
generates  a  sound  wave  of  wavelength  A.  Within  a  distance 
a  /A  from  the  emitter,  Fresnel  diffraction  takes  place, 
and  the  sound  beam  is  cylindrical  with  a  radius  a.  Within 
that  distance,  the  amplitude  of  the  sound  beam  averaged  over 
a  cross-section  of  the  cylinder  decreases  as  e  ,  where 
a  is  the  absorption  coefficient  of  the  medium,  and  x  is 
the  distance  from  the  source.  Standing  waves  are  a  problem, 
although  at  high  enough  absorption  and  at  large  distances, 
their  effect  becomes  negligible.  Mechanical  and  electrical 
cross-talks  also  give  difficulties  if  they  are  not  neutralized 
properly.  In  gases,  this  method  has  been  applied  quite 
successfully  at  3  Mc/sec/atm  and  higher,  in  particular, 

r-7  Q 

by  Meyer  and  Sessler  ,  and  by  Sluijter,  Knaap  and  Beenakker. 

All  methods  making  use  of  continuous  waves  present 
the  following  disadvantages:  standing  waves  can  affect  the 
results,  and  moreover,  if  power  is  continuously  dissipated 
in  the  form  of  heat  in  the  medium,  the  temperature  of  the 
fluid  may  be  altered  locally,  and  measurements  taken  at  a 
temperature  higher  than  that  recorded.  In  the  pulse  method. 
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these  problems  do  not  arise:  if  the  pulse  repetition  rate 

is  low  enough,  standing  waves  cannot  be  set  up,  and  if  the 

duty  ratio  is  of  the  order  of  1:500,  the  heating  effect  is 

reduced  in  the  same  proportion.  The  pulse  method  also  permits 

to  distinguish  in  time,  between  the  received  pulse  and  other 

unwanted  signals.  One  objection  may  be  raised  against  the 

pulse  method  however;  since  a  short  pulse  contains  a  whole 

spectrum  of  frequencies,  and  since  the  sound  absorption 

varies  with  the  frequency,  there  is  an  error  involved  in 

q 

taking  measurements  with  short . pulses «  Pellam  and  Galt 
demonstrated  that  when  the  absorption  coefficient  varies  as 
the  square  of  the  frequency,  and  when  the  pulse  contains 
several  oscillations,  the  error  involved  is  equal  to  the 
square  of  the  reciprocal  of  the  number  of  oscillations  per 
pulse.  For  example,  when  a  pulse  contains  ten  oscillations, 
the  error  is  one  part  in  a  hundred,  which  is  usually  less 
than  the  errors  associated  with  the  measurement  of  other 
parameters,  such  as  the  pressure,  the  voltage  etc.  The  main 
disadvantage  of  the  pulse  method  is  the  poor  signal -to-noise 
ratio  associated  with  the  large  bandwidth  of  the  receiver. 

In  conclusion  ,  no  method  is  obviously  the  best  one. 
A  compromise  must  be  reached,  which  depends  on  the  problem 
to  be  solved,  as  well  as  on  the  facilities  available  in  a 
particular  laboratory. 
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2.  Apparatus 

In  order  to  observe  a  single  relaxation  process  in 
hydrogen,  measurements  have  to  be  taken  at  temperatures  below 
150°K.  Since  liquid  nitrogen  was  the  most  readily  available 
refrigerant,  a  gas  vessel  was  constructed  which  is  suitable 
for  that  temperature.  Another  requirement  imposed  on  the 
apparatus  was  the  range  of  f requency-to-pressure  ratios  (f/p) 
at  which  experiments  were  to  be  conducted.  Since  the  maximum 
of  the  absorption  peak  for  hydrogen  occurs  at  around  10  Mc/sec/atm, 
the  apparatus  was  built  to  cover  an  appropriate  range  of  f/p: 
from  3  to  50  Mc/sec/atm.  Equations  (1-80),  (l-8l)  and  (1-82) 
show  that  t  always  appears  in  the  combination  .  Assuming 

that  rotational  energy  transfer  is  due  to  binary  collisions 
only,  f  is  inversely  proportional  to  the  collision  rate,  and 
consequently  inversely  proportional  to  the  pressure.  The 
absorption  can  thus  be  measured  as  a  function  of  f/p,  rather 
than  of  The  proportionality  factor  between  f/p  and 

is  independent  of  frequency  and  pressure.  Since  it  is  easier 
experimentally  to  vary  the  pressure  rather  than  the  frequency, 
measurements  were  taken  at  a  fixed  frequency,  for  various 


pressure  values. 
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A.  Mechanical  part. 

A  schematic  diagram  of  the  cryostat  appears  in 
fig.  4.  The  cryostat  is  supported  vertically  and  can  thus  be 
immersed  in  a  refrigerant  bath.  The  double  bellows  arrangement , 
adapted  from  Stewart,10  has  two  advantages:  firstly,  the 
volume  remains  constant  while  the  distance  between  the  two 
transducers  is  varied,  and  secondly,  this  arrangement  consti¬ 
tutes  a  temperature-resistant  motion  seal  --  other  types  of 
motion  seals  while  efficient  at  room  temperature,  may  not  be 
leak-proof  at  liquid  nitrogen  temperature.  The  elasticity  of 
the  bellows  permits  the  length  of  the  sound  path  to  be  varied 
from  1.5  to  14  mm.  The  position  of  the  movable  transducer  is 
determined  by  means  of  a  micrometer. 

Indium  O-rings  are  used  as  seals  between  the 
stainless  steel  flanges  which  are  soldered  to  the  bellows, 
and  the  flanges  of  the  central  rigid  portion  of  the  cryostat. 

Spring  loaded  adjustment  screws  permit  the  transducers 
to  be  set  as  parallel  as  possible.  The  right  degree  of  parallel¬ 
ism  is  achieved  when  the  received  signal  is  maximum.  It  is 
also  important  that  this  parallelism  be  maintained  as  the 
movable  transducer  is  brought  forward  or  backward:  This  is  the 
purpose  of  the  driving  shaft  fixed  to  the  upper  plate  of  the 
cryostat,  which  slides  inside  a  tight-fitting  bushing. 


■ 
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Pig.  4 

Cross  Section  of  Apparatus 


Micrometer 


Supporting 

Plate 


Sliding 

Shaft 


Guiding 

Bushing 


BNC 

Connector 


Transducer 

Assembly 


Indium 

O-ring 


Gas  Line 


Window 


Transducer 

Assembly 


Indium 
0— ring 


Levelling 

Screws 


Electrical 

Wire 


Brass 


Stainless 

Steel 


51 


B.  Mounting  of  the  Transducers 

The  transducers  are  X-cut  quartz  crystals  *  one  inch 
in  diameter,  with  a  fundamental  frequency  of  3.00  Mc/sec. 

They  are  cemented  to  a  block  of  Lucite.  A  cross-section  of 
the  transducers  support  appears  in  fig.  5.  Many  attempts 
were  made  at  providing  the  transducers  with  proper  electrodes. 
Films  of  copper,  aluminum  and  chromium  were  deposited  onto 
the  surfaces  of  the  transducers  by  evaporation  in  vacuo.  The 
electrodes  thus  obtained  were  too  thin  and  did  not  adhere 
satisfactorily  to  the  quartz  substrate.  Coating  the  trans¬ 
ducers  with  silver  metallizing  paint  (electrically  conducting), 
proved  to  be  successful  to  some  extent,  although  the  paint 
had  a  tendency  to  dissolve  in  the  grease  used  to  bond  the 
transducers  to  the  Lucite.  Another  attempt  consisted  of 
cementing  thin  discs  of  aluminum  foil  onto  the  surfaces  of  the 
transducers.  This  also  was  a  rather  satisfactory  arrangement, 
although  the  transducers  seemed  to  be  too  heavily  loaded,  with 
the  result  that  their  resonant  frequency  was  lowered.  Conse¬ 
quently,  the  following  compromise  was  reached:  the  surface 
of  the  transducer  which  is  in  contact  with  the  gas  is  coated 
with  silver  paint,  while  the  surface  which  is  cemented  to  the 
Lucite  is  covered  with  an  electrode  cut  out  from  a  sheet  of 
thin  aluminum  foil.  Silicone  vacuum  grease  proved  to  be  the 
best  bonding  agent.  At  liquid  nitrogen  temperature,  the  resonant 
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Fig.  5 

Transducers  Mount 
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frequency  of  the  transducers  is  2.85  Mc/sec.  The  lowering  of 
the  resonant  frequency  is  seemingly  due  to  the  extra  thickness 
of  the  foil  and  the  thin  film  of  grease  used  to  bond  the 
electrode  to  the  quartz  crystal. 

Unloaded  quartz  crystals  have  a  very  high  "Q"  factor, 
which  can  reach  10,000.  Such  a  highly  tuned  device  cannot 
respond  to  a  pulse  of  short  duration.  In  the  pulse  method, 
it  is  thus  necessary  to  reduce  the  "Q"  factor  of  the  transducer, 
or  in  other  words,  to  increase  its  bandwidth.  This  is 
achieved  by  backing  the  transducer  with  a  high  impedance  and 
high  loss  material.  Lucite  satisfies  these  two  conditions, 
and  is  used  to  load  the  transducers.  Cemented  to  a  block  of 
Lucite,  the  quartz  crystals  have  a  "Q"  of  approximately  15. 

Electromagnetic  shielding  of  transducers  calls  for 
coaxial  types  of  connectors.  It  is  for  this  reason  that  the 
arrangement  shown  in  fig.  5  was  adopted:  the  ground  connection 
is  provided  by  a  brass  hollow  cylinder  surrounding  the  cylin¬ 
drical  block  of  Lucite  through  the  centre  of  which  passes  the 
"live"  connector  -  a  thin  brass  rod  terminated  in  a  spring 
loaded  plug.  One  end  of  the  brass  cylinder  is  such  that  it 
covers  the  periphery  of  the  transducer’s  face,  thus  providing 
electrical  contact  to  the  ground. 
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C.  Vacuum  system 

The  gas  vessel  is  evacuated  by  means  of  a  mechanical 
pump  capable  of  an  ultimate  vacuum  of  10  mm  Hg.  As  a 
preliminary  measure,  the  apparatus  was  checked  for  leaks,  and 
connected  to  a  mercury  diffusion  pump.  After  evacuation  by 
means  of  the  diffusion  pump,  the  vessel  maintains  overnight  a 
vacuum  better  than  10  mm  Hg. 

Pressures  are  read  from  a  mercury  manometer  with 
an  accuracy  of  -  0.5  mm. 

D.  Electronic  Circuitry 

A  block  diagram  of  the  electronic  components 
appears  in  fig.  6. 

A  R.F.  pulse  of  adjustable  duration,  frequency 
and  repetition  rate  originates  from  an  oscillator  driven  by 
a  blocking  oscillator.  Matching  of  the  transducer  to  the  output 
of  the  oscillator  is  effected  by  means  of  a  tuned  air-core 
transformer  with  variable  number  of  turns  and  variable  turns  ratio. 

The  receiving  transducer  is  connected  to  a  preamplifier 
with  tunable  input  coils.  The  signal  is  then  amplified  in 
a  wide -band  amplifier  with  a  maximum  gain  of  80  db.  The  video 
output  signal  from  the  amplifier  is  measured  by  means  of  a 
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Fig.  6 


Electronic  Circuitry 
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"sampling"  circuit  adapted  from  Kamm  and  Bobm.11  The  output 
of  the  sampling  circuit  is  recorded  on  a  2.5  mv  full-scale 
detector . 

Such  a  method  requires  that  the  gain  of  the  amplifier 
be  very  stable;  for  this,  all  filaments  of  the  amplifying 
system  were  heated  by  large  capacity  batteries  connected  to  a 
battery  ~  eliminator. 

3.  Procedure 

The  cryostat  is  filled  with  the  gas  to  be  examined,, 
to  a  pressure  slightly  above  atmospheric  pressure.  The 
received  signal  is  recorded  at  various  positions  of  the 
movable  transducer.  A  typical  record  appears  in  fig.  J .  Each 
record  contains  in  general  fifteen  double  peaks,  corresponding 
to  fifteen  different  transducers  separation.  The  peaks  are 
double  since  each  pulse  is  "sampled"  twice,  backward  and 
forward.  The  amplitudes  are  then  read  from  these  records  and 
the  natural  logarithm  of  these  amplitudes  are  plotted  versus 
the  distance.  The  absorption  coefficient  is  determined  from  the 
slope  of  these  straight  lines.  A  typical  graph  of  In  A  versus 
x  is  shown  in  fig.  8. 

The  pressure  is  then  decreased  to  the  next  lower  value, 
and  the  above  procedure  is  repeated,  similarly  down  to  the 
lowest  workable  pressure. 


■ 
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Fig.  7 


Typical  Record 


\ 
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Fig .  8 


Curve  of  In  A  vs  x 


CHAPTER  III 


RESULTS  AND  DISCUSSION 


1.  Sound  Absorption  Measurements  in  Helium  at  Room  Temperature 

Before  undertaking  measurements  in  hydrogen,  it  was 
necessary  to  test  the  apparatus.  For  this,  measurements  were 
first  taken  in  helium  at  room  temperature.  The  results  appear 
in  fig.  9  and  table  2.  They  are  compared  with  the  theoretical 
curve  calculated  from  equation  (I-58) 

.  -II 

<A  class  -  g  R 
2 

where  R  =  Po  V o  ,  and  V  =  / VpQ 

t[  u> 

The  data  necessary  for  calculating  ^qass  appear 
in  table  1.  The  value  of  the  viscosity  was  calculated  from 
Keesom's  formula 

tl  =  5.023  T0,647 

in  which  the  units  are  in  micropoises  if  T  is  in  degrees 
Kelvin.  The  results  are  in  good  agreement  with  the  theoretical 
values.  The  estimated  error  is  certainly  less  than  10$  even 
at  high  f/p  values. 

2.  Sound  Absorption  Measurements  in  Helium  at  J6.J  °K. 

The  apparatus  was  next  tested  for  proper  operation 
at  liquid  nitrogen  temperature.  Due  to  the  low  atmospheric 
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Table  2 

Absorption  in  Helium  at  296°K. 
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Table  2 


1/p 

f/p 

o( 

<*/f 

Atm-1  " 

Me/  / 

cm"'*' 

( cm  Me 

/sec/atm 

1.05 

3.15 

0.48 

0.156 

1.34 

4.02 

0.63 

0.210 

1.53 

4.59 

0.805 

0.264 

1.96 

5.88 

0.97 

0.324 

2.57 

7.71 

X) 

OJ 

• 

1 — 1 

0.430 

3.08 

9.24 

1.52 

0.510 

3.70 

i — 1 

• 

i — I 
i — 1 
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Figure  9 

<*/f  vs  f/p  for  helium  at  296  °K 
0  Experimental  points 


Theoretical  curve 
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pressure  in  this  laboratory,  liquid  nitrogen  boils  at  76.7°  K. 

The  results  in  helium  at  liquid  nitrogen  temperature 

appear  in  fig.  10  and  the  individual  values  figure  in  table  3. 

Again,  they  are  compared  with  the  theoretical  values.  The  value  of 

the  viscosity  was  taken  this  time  from  Corremans,  van  Itterbeek, 

12 

Beenakker,  Knaap  and  Zandbergen  .  See  table  1.  The  general  accuracy 
is  again  better  than  10$. 


3.  Sound  Absorption  Measurements  in  Normal  Hydrogen  at  76.7°  K. 


Although  both  para  and  ortho  modifications  are  present  in 
normal  hydrogen,  the  gas  should  exhibit  a  single  relaxation  process 
at  liquid  nitrogen  temperature,  since  the  first  excited  level  of 
ortho  hydrogen  is  not  populated  at  that  temperature.  As  can  be 
seen  from  fig.  11,  and  table  4,  the  results  obtained  in  this  labo¬ 
ratory  indeed  fit  a  single  relaxation  curve,  within  the  experimental 
accuracy. 

The  theoretical  curve  was  calculated  from  the  following 
expression 


4/t  _  1  ( !-€.  u££X—J\ 

'tI  nsj 


II  £  uJ  L 


k\_ 


l  +cu  ti 


PS 

which  is  obtained  from  equations  (l-8o)  and  (l-8l).  This  curve  has 


at 


C"  a  ~  ( 1  (Herzfeld  and  Litovitz)1^ 


pS 


The  relaxation  strength  was  calculated  from 


£ 


c'  (<h°  -c°v) 


a  maximum 


a  Isnraofl  nl 
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Table  3 

Absorption  in  helium  at  76.7°K. 


67 


Table  3 


1/p 
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Figure  10 


vs  f/p  for  helium  at  76.7 
Experimental  points 


Theoretical  curve 
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Table  4 

Absorption  in  normal  hydrogen  at  76.7°  K 
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Table  4 


f/p 

class 

£^exp 

y 

<d 

d/f 

V..e/at. 

-1 

cm 

-1 

cm 

-1 

cm 

( cm  Me 

2.42 
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Figure  11 

c</ f  vs  f/p  for  normal  hydrogen  at  76.7°  K 

©  Experimental  points 

-  -  -  -  Theoretical  curve  calculated  for  £  =  0.02 
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The  values  of  the  specific  heats  were  taken  from 

14  0 

Farkas  ,  See  table  1.  For  normal  hydrogen  at  76.7  6  =  0.02. 

The  value  of  °(^was  obtained  by  subtracting  &  ,  from 

ClaSS 

the  experimentally  determined  value  of  c\  which  is  denoted  by  cK 

exp 

Since  hydrogen  is  not  a  monatomic  gas.,  equation  (l-59)  rather 
than  (1-58)  had  to  be  used  in  calculating  : 

C  _Lcl  S  S 


The  data  for  the  viscosity  were  taken  from  Coremans  et  al., 

while  the  data  for  the  thermal  conductivity  were  taken  from 
15 


Ubbink' 


See  table  1. 


In  general,,  the  accuracy  of  the  results  seems  poorer 

than  in  the  case  of  helium.  It  is  reasonable  however,  to  assume 

that  was  determined  to  the  same  degree  of  accuracy  as  in  the 

case  of  helium.  The  higher  discrepancy  between  the  theoretical 

'  /  ‘ 

and  the  experimental  values  of  cK/f  is  due  to  the  fact  that 
represents  only  a  fraction  of  the  total  absorption.  Therefore, 
anjerror  of  10$  on  constitutes  a  higher  percentage  error  on 

o{  .  For  example,  at  20.3  Mc/sec-atm  ,  o{  was  found  equal  to 

0X  p 

1.5^.  A  10$  error  is  0.15.  Since  o<;  is  obtained  by  subtracting 


from  ^  is  0.8l.  The  absolute  error  0.15  represents 

a  relative  error  of  more  than  18$  on  .  As  the  classical 


°^class 


absorption  becomes  more  important, cK  represents  a  smaller  fraction 

> 

of  the  total  absorption  and  the  accuracy  to  which  <A  can  be  deter¬ 


mined  becomes  lesser.  Taking  this  fact  into  account,  the  results 


.  t  . 

•50 

-u  **d 
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in  fig. 11  agree  with  the  theoretical  curve  well  within  the 
experimental  error. 

The  maximum  of  the  theoretical  curve  which  best  fits  the 
experimental  points  is  located  at  7.65  Me  /sec-atm.  At  this  point, 
OU^pS  »  (l  -  0.02)-1///\  Since  is  the  adiabatic  relaxation  time 

at  unit  pressure,  the  maximum  of  the  relaxational  absorption 
occurs  at  a  frequency 

27 r  X  7.65  X  106  sec'1. 


From  this,  f  g  is  found  to  be  equal  to 


(1  -  0.02) 


-1/4 


2.11  X  10' 


•  8 


sec . 


2  7f  X  7.65  X  10 


6 


The  isothermal  or  true  relaxation  time  can  be  calculated  from  the 
following  relation 


Tps 

t  pT 

With  C°  =  2.58R 
P 


C 


OO 


c 


£ 

o 

P 


OO 


and  C  r  2.50R, 

p 


f-  _Q 

l  T  =  2.17  x  10  see 


4.  Sound  Absorption  Measurements  In.  Para-enriched  Hydrogen. 

An  attempt  was  made  at  producing  pure  para  hydrogen. 

Normal  hydrogen  was  liquefied  in  the  presence  of  a  catalyst  (hydrous 
ferric  oxide) .  Since  liquid  hydrogen  is  not  available  in  this 


.  V'-/'.  ,  •  •  •• 


■ 


.  ■ 
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laboratory,  conversion  was  effected  at  a  temperature  in  the  neighbour¬ 
hood  of  20°  K,  above  the  surface  of  the  liquid  inside  a  ewar 
containing  liquid  helium.  A  few  reasons  could  have  contributed  to 
the  lack  of  success  of  this  experiment;  the  catalyst  for  instance 
was  probably  not  activated  sufficiently,,  and  the  liquified  hydrogen 
may  not  have  been  in  contact  with  the  catalyst  for  a  long  enough 
time . 

Nevertheless,  a  definite  enrichement  in  para  content  was 
achieved,  and  sound  absorption  measurements  were  taken  in  this  para- 
ortho  mixture.  The  results  are  shown  in  fig.  4  and  the  individual 
values  figure  in  table  5. 

)  . 

From  the  maximum  value  of  o(/f,  the  relaxation  strength 
was  found  to  be  equal  to  0.045.  The  theoretical  values  of  the 
relaxation  strength  of  normal  and  pure  para  hydrogen  are  0.02  and 
0 .06  respectively.  A  relaxation  strength  of  0.045  corresponds  to 
a  rotational  specific  heat  C1  of  0.19  R. 


The  maximum  of  the  absorption  curve  occurs  at  a  f/p  value 


of  8. 

1  Mc/sec-atm 

.  From 

this,  ^  .  is  found  to  be 

equal  to 

1.96 

g 

X  10  sec. 

The 

isothermal  relaxation  time 

in  this  case  is 

Tpt  - 

1.96 

X 

0 

2.69  X  10-  sec 

2.50 

or 

h3 

II 

2.11 

X 

n  —  8 

10  sec. 

which  agrees  within  a  few  percent  with  the  isothermal  relaxation 

time  for  the  0  - >  2  transition  found  in  normal  hydrogen. 

It  is  of  interest  to  compare  the  present  results  with 
those  obtained  by  other  authors.  Sound  absorption  measurements 
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Table  5 

Absorption  in  para-enriched  hydrogen  at  76.7°  K 
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Table  5 


f/p 

^class 

^exp 

j 

0< 

o(J/f 

Me 

sec  atm 

cm-^ 

-1 

cm 

-1 

cm 

( cm  Me 

3.03 

0.109 

1.72 

1.6l 

0.57 

4.28 

0.154 

2 . 40 

2.25 

0.79 

5.90 

0.214 

2.68 

2.47 

0.87 

7.10 

0.258 

2.93 

2.67 

0.94 

8.10 

0.293 

3.18 

2.89 

1.01 

9 . 60 

0.345 

3.08 

2.74 

0.96 

10.6 

0.375 

3.08 

2.71 

0.95 

11.5 

0.415 

3.10 

2.69 

0.92 

12.9 

0.465 

2.83 

2.37 

0.83 

14.6 

0.530 

2.98 

2.45 

0 .86 

16.9 

0.610 

2.84 

2.23 

0.78 

24.6 

0.892 

2.67 

1.78 

0.62 

30.9 

1 — 1 

1 — 1 

• 

1 — 1 

2.33 

OJ 

OJ 

• 

1 — 1 
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Figure  12 

d/f  vs  f/p  for  para-enriched  hydrogen  at  76.7°  K 
0  Experimental  points 

-  -  -  -  Theoretical  curve  calculated  for  £-0.045 


78 


in  hydrogen  have  been  taken  by  several  authors  for  many  years. 

In  1937  and  1938.,  van  Itterbeek  and  Mariens1^,  and  van  Itter- 

17 

beek  and  Thys  reported  results  obtained  with  an  ultrasonic 

interferometer  in  the  region  where  is  constant, 

class 

i.e.  at  f/p  values  well  below  the  region  where  ^ 1 /f  goes 
through  a  maximum.  These  results  however  did  not  show  the  right 
variation  with  pressure. 

l8 

Stewart  and  Stewart1  conducted  measurements  in  normal 

hydrogen  with  an  ultrasonic  interferometer  at  a  frequency  of 

1.95  Mc/sec.  With  an  accuracy  of  20$,  their  results  indicated 

that  the  maximum  of  the  absorption  curve  occurred  at  a  f/p 

value  of  the  order  of  10  --  12  Mc/sec/atm. 

Improved  measurements  taken  by  van  Itterbeek  and 

Verhaegen1^  at  79.5  °K  and  90  °K  showed  the  right  variation 

with  pressure.  These  results  obtained  again  in  the  region  where 

<^/cK  i  is  constant,  yield  a  value  of  approximately 

'  c  -La  s  s 

-8 

2.5  X  10  sec.  It  can  be  shown  that  in  that  region, 

class  =  1  +  0.067_tpT 

to 

where  t  is  the  average  time  between  collisions  at  a  pressure 
of  one  atmosphere. 

20 

Using  the  pulse  method,  Parbrook  and  Tempest  took 
measurements  in  normal  hydrogem  at  298  °K.  The  results  show  a 
lower  and  broader  curve  which  does  not  fit  a  single-relaxation 
curve.  This  is  to  be  expected  since  at  room  temperature, 
normal  hydrogen  has  five  excited  levels  giving  rise  to  four 


sarfXo 

■ 

.  jo 
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possible  transitions,  with  possibly  four  different  relaxation 
times . 

All  measurements  made  in  either  pure  para  hydrogen 
or  in  normal  hydrogen  above  150  °K  have  to  be  interpreted  in 
terms  of  two  or  more  relaxation  processes.  This  necessitates 

the  use  of  many  assumptions.  In  order  to  analyse  the  results 

o  21 

which  he  obtained  at  293  K  in  normal  hydrogen,  Sessler  made 

the  following  assumptions:  the  relaxation  time  for  the  1  - ->  3 

transition  is  1.5  times  that  for  the  0  - ->2  transition,  and 

the  relaxation  time  for  the  3  - ^5  transition  is  the  same 

as  that  for  the  1  - ->  3  transition.  He  assumed  moreover 

that  all  transitions  occur  independently.  From  this,  Sessler 

found  that  at  room  temperature  the  isothermal  relaxation  time 

for  the  0  - ->2  transition  (fQ2)  is  1.29  X  10-^  sec. 

22 

Geide  who  took  measurements  in  pure  para  hydrogen 

at  293  °K,  also  assumed  that  the  2  - >  4  transition  occurs 

independently  of  the  0  - }  2  transition.  With  this  assumption, 

his  measurements  yield  a  value  of  Tq£  =  1*21  X  10-^  sec. 

It  must  be  emphasized  that  these  relaxation  times 

-8 

(1.29  and  1.21  X  10  sec  )  were  obtained  from  measurements 
taken  at  room  temperature,  and  it  is  reasonable  to  expect  that 
they  should  differ  from  the  relaxation  times  obtained  from 
measurements  taken  at  lower  temperatures. 

The  following  considerations  will  give  an  idea  of  the 
way  in  which  the  relaxation  times  can  be  expected  to  vary 
with  temperature. 
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As  was  pointed  out  at  the  end  of  the  first  chapter, 
the  isothermal  relaxation  time  at  unit  concentration  is  the 
reciprocal  of  the  sum  of  the  reaction  rates.  In  para  hydrogen 


-1 

T 


r02  * 


R, 


20 


The  relation  between  Rq2  and  R2Q  is  provided  by  the  condition 
of  statistical  equilibrium, 

,g  ,  E/«T 

if  • 


R 


20  = 


R 


02 


where  gQ  and  g2  are  the  statistical  weights  of  the  rotational 

levels  J  =  0  and  J  =  2  respectively.  E  is  the  energy  difference 

between  the  energy  of  the  state  J=2  and  J  =  0. 

Now,  the  reaction  rate  is  proportional  to  the  number 

of  molecules  striking  unit  area  per  unit  time  at  pressure  p  and 

temperature  T.  From  kinetic  theory,  this  number  of  molecules 

is  /  ^  }  1/2  p.  The  reaction  rate  is  also  proportional 

\fTmkT  / 

to  the  efficiency  of  the  collisions  in  inducing  rotational 
energy  transfer.  The  efficiency  of  the  collisions  is  expressed 

in  terms  of  an  averaged  inelastic  cross  section  for  the  0  - ->2 

transition  at  temperature  T  and  pressure  p.  This  cross  section 
is  denoted  by  Qq2/  t.  Therefore, 


R 


8 


1/2 


02 


T T  mkT 


p  <  q02)  t  , 


and 


.  r-fOO  P^>  **  '  " 


. 
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^  ^02^  T  is  stained  Sy  averaging  Q(K),  the  Inelastic  cross 

section  corresponding  to  a  wave  number  K,  over  the  Maxwellian 

distribution  of  wave  numbers  K  at  temperature  T. 

2  c  26  27 

Various  authors  ^  *  have  calculated  Q(K)  for 

different  forms  of  the  intermolecular  potential.  Using  a 

25 

Morse-type  potential  with  two  angle-dependent  terms.,  Davison 
has  computed  the  values  of  Q(K)  as  function  of  K  .  His  results 
may  be  approximated  as  follows: 


Q(K)=  0  (K2  -K2  ) 


where  c  =  0.043  and  =  6  in  atomic  units.  In  this  system 
of  units j  Planck’s  constant  is  equal  to  2T T. 

Averaged  over  the  Maxwellian  distribution  of  wave 
numbers  at  temperature  T,  Q(K)  yields 


<Q02>  t  =  exp  (  -E/kT  )  4cmkT  (  1 


Inserting  < Q  ? V  into  the  expression  one 


02/T 


gets 


l  ~ 1 

ip  — 


8 


1/2 


JT  mkT 


o  exp  (-  E/kT  ) 


4cmkT  (  1  + 


Above  150  °K*  the  factor  T  * 


is  dominant  and  one 


should  expect  the  relaxation  times  to  decrease  as  the  temperature 
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increases . 

A  rough  calculation  shows  that  the  relaxation  time 


at  room  temperature  should  be  approximately  half  the  relaxation 
time  at  liquid  nitrogen  temperature. 


Indeed, 


at  80  °K,  E/kT  1 

and 

e'1  ^  0.37 

at  300  °K,  E/kT  -  0.3 

and 

e~0,3  ~  0.72 

g0-  1  and  g2=  3. 

Therefore, 


Co.  33  f  0.37)  x  (if  0.5) 
(o.33  +  0.72)  (l  to. 15) 


1 300 

r  80 


1/2 


The  relaxation  time  deduced  from  Sessler's  measurements 

-8 

at  room  temperature  (1.2  X  10  sec)  is  roughly  half  the 

relaxation  time  deduced  from  measurements  at  liquid  nitrogen 

-8 

temperature  (2.1  X  10  sec). 

The  computation  of  /  Qq2")  ^  involves  many  assumptions 
and  especially,  the  use  of  certain  adjustable  parameters.  One 
may  question  the  validity  of  the  assumptions  made,  as  well  as 
the  values  of  the  parameters  used.  The  best  way  of  verifying 
the  reliability  of  the  theory  consists  of  comparing  the  theore¬ 
tical  and  the  experimental  results  over  a  wide  range  of  tempe¬ 
rature.  In  that  respect,  is  not  the  best  quantity  to/konsider 

for  making  comparisons,  because  the  temperature  dependence  of 
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'f¬ 
lu  T  is  affected  by  the  fact  that  i T  depends  on  the  number  of 

molecules  hitting  unit  area  per  unit  time.  As  seen  previously, 
the  relaxation  time  is  inversely  proportional  to  the  number 
of  molecules  hitting  unit  area  per  unit  time;  it  is  thus  propor¬ 
tional  to  the  mean  time  between  collisions  T  .  At  pressure  p, 

c 

te1  =  ™-2  (JLA  1/2  p 

V7rmkT  / 


where  <r  is  the  diameter  of  the  molecules  treated  as  hard  spheres. 

Consequently,  the  collision  number  Z  =  is 

better  suited  for  comparisons  between  theory  and  experiment. 


Using  the  expressions  for  t 


-1 

T 


,-l 


IT  (Tz 


1  + 


and  t “1  , 
g  E/kT 


Since  <(  Qq2)  rp  results  from  averaging  Q(K)  over  the  Maxwellian 
distribution  of  the  wave  numbers  at  temperature  T,  it  can  be 
written  as 


<Q02>  T  =  exp  (-  E/kT)  l(T) 


where  l(T)  stands  for  an  integral  which  depends  on  the  form 
of  the  intermolecular  potential.  With  <(  Qq2)  ip  expressed 
as  above. 


Z 


-1 


1 

Tr<r  2 


g 

■zr-  -H  exp  (-E/kT) 
S2 


I  ( T) 


Therefore,  the  temperature  dependence  of  Z  is  seen 
to  be  mostly  determined  by  the  temperature  dependence  of  l(T) . 


■ 
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Taking  l(T)  —  4cmkT  (l  +■  E/OkT^  as  before. 


(0.33  +  0.72) 

-  X 

(0.33  +  0.37) 


(1+0.15)  300 

-  X  - 

(1+0.5  )  80 


Z300 


This  value  is  in  good  agreement  with  the  experimental 

results : 


at 

300  °K, 

t  = 

0  c 

0.7 

X 

10  10  sec. 

at 

00 

0 

0 

t  - 

n 

0.3 

X 

in"10 

10  sec. 

Consequently,  using  Sessler's  value  of  the  relaxation 
time  at  room  temperature. 


Z 


300 


1.2  X 

0.7  X 


170  collisions. 


At  liquid  nitrogen  temperature, 

_  o 

pi  y  in  0 

Z0n  ^  —  700  collisions. 

0.3  x  10  u 

These  values  of  the  collision  numbers  are  in  good 

agreement  within  an  order  of  magnitude,  with  the  theoretical 

24 

values  deduced  from  the  calculations  of  Sluijter  et  al.  ,  and 

30 

those  of  Brout 

The  only  results  obtained  in  hydrogen  below  150  °K, 

which  are  reported  in  the  literature,  are  those  of  Sluijter, 

8  23  24 

Knaap  and  Beenakker  *  ’  .  From  their  measurements  taken  in 
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pure  para  hydrogen  at  77.3  °K  and  90  °K,  the  values  of  £ 

-ft  -ft 

are  found  to  be  2.2  X  10*"  sec  and  2.1  X  10-  sec 
respectively. 

Measurements  by  Sluijter  et  al.  In  normal  hydrogen 
at  77.3  °K,  however  indicated  that  the  maximum  of  the  relaxa- 
tional  absorption  occurred  at  approximately  11  Mc/sec/atm  or 
even  higher.  This  indicated  a  value  of  £"  =■  1.^  X  10~^  sec. 

This  important  difference  between  the  relaxation  time 

of  para  hydrogen  at  100$  concentration,  and  the  relaxation  time 

-ft  ft 

at  25$  concentration  (2.2  X  10~  sec  and  1.4  X  10”  sec 

respectively.  ),  raises  the  question  as  to  whether  £"02  is  the 

same  when  collisions  occur  between  para  molecules  only,  as  it 

is  when  para  and  ortho  molecules  collide  together.  A  shorter 

relaxation  time  in  the  case  of  normal  hydrogen  would  mean  that 

ortho  molecules  are  more  efficient  than  para  molecules  at 

inducing  rotational  energy  transfer  during  a  collision. 

The  present  results  seem  to  indicate  that  such  is  not 

A- 

the  case,  since  the  values  of  obtained  from  measurements 

in  normal  hydrogen  are  in  good  agreement  with  the  value  of  c02 
obtained  by  Sluijter  et  al.  from  measurements  in  pure 
para  hydrogen. 


■ 

i: 


■ 

. 
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CONCLUSION 


By  means  of  the  pulse  technique  at  liquid  nitrogen 
temperature,  the  sound  absorption  was  measured  in  normal  and 
para-enriched  hydrogen  at  frequency-to-pressure  ratios  in  the 
range  3  to  50  Mc/sec/atm. 

These  measurements  are  in  good  agreement  with  the 
theory  first  proposed  by  Herzfeld  and  Rice.,  both  as  to  their 
dependence  on  f/p  and  their  absolute  magnitude. 

The  peak  in  the  absorption  curve,  referred  to  a 

pressure  of  one  atmosphere,  occurs  when(JU-1=f'  ,  the  adiabatic 

po 

relaxation  time.  From  this,  =Tq2  was  deduced,  and  the 

'N  -8 

value  obtained  was  Lq 2=2.1  X  10  sec,  for  both  para-enriched 
and  normal  hydrogen. 

This  value  of  C"q2  is  an  good  agreement  with  that 
deduced  from  measurements  taken  by  the  direct  method  at  the 
same  temperature,  in  pure  para  hydrogen.  It  differs  however 
from  the  value  of  i  deduced  from  measurements  taken  by 
other  authors  in  normal  hydrogen  at  the  same  temperature. 

Although  no  claim  is  made  at  a  particularly  high 
experimental  accuracy,  the  fact  that  the  number  of  measurements 
(i.e.  the  number  of  points  on  the  absorption  curve)  taken  in 
normal  hydrogen  at  liquid  nitrogen  temperature  and  reported 
in  this  thesis,  is  considerably  larger  than  the  number  of 
similar  measurements  taken  and  reported  by  any  other  authors. 


may  confer  to  the  present  results  a  certain  degree  of  reliability. 
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The  relaxation  theory  of  Herzfeld  and  Rice  is  a 
thermodynamical  one  and  as  such  can  only  provide  a  macroscopic 
description  of  the  phenomenon  involved.  A  comparison  of  the 
relaxation  times  obtained  at  liquid  nitrogen  temperature,  with 
those  obtained  by  other  authors  at  room  temperature  shows  that 
the  microscopic  theory  predicts  quite  satisfactorily  the  tempera¬ 
ture  variation  of  the  relaxation  times.  A  full  verification 
of  the  microscopic  theory  must  however  await  experimental 
results  at  other  temperatures.  Moreover,  theories  are  not 
precise  enough  at  the  present  time  to  give  absolute  magnitude 
of  the  relaxation  times  to  better  than  within  an  order  of 
magnitude . 

As  yet,  no  reliable  information  has  been  obtained 

about  the  relaxation  time  of  the  1  - }  3  transition  in  ortho 

hydrogen.  If  measurements  could  be  performed  in  pure  ortho 

hydrogen,  the  results  would  be  most  interesting.  Unfortunately, 

29 

the  production  of  ortho  hydrogen  although  feasable  9  remains 
difficult  especially  in  view  of  the  large  amount  of  gas  neces¬ 
sary  to  fill  a  cryostat.  Measurements  in  normal  hydrogen  or 
in  various  mixtures  of  ortho  and  para  hydrogen  between  140  °K 
and  l80  °K  could  also  provide  some  information  about  the  relaxa¬ 
tion  time  of  the  ortho  modification.  Results  in  ortho  hydro¬ 
gen  would  constitute  another  way  of  verifying  the  microscopic 


theories . 


■  -  ■  f 

. 
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